fitting polynomial regression with r